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Abstract. We present a unified general method for the asymptotic study of 
graphs from the so-called "subcritical" graph classes, which include the classes 
of cacti graphs, outerplanar graphs, and series-parallel graphs. This general 
method works both in the labelled and unlabelled framework. The main results 
concern the asymptotic enumeration and the limit laws of properties of random 
graphs chosen from subcritical classes. We show that the number gn/n] (resp. 
gn) of labelled (resp. unlabelled) graphs on n vertices from a subcritical graph 
class Q = UnSn satisfies asymptotically the universal behaviour 

g„ = cn-^/S"(l + o{l)) 

for computable constants c, 7, e.g. 7 9.38527 for unlabelled series-parallel 
graphs, and that the number of vertices of degree k {k fixed) in a graph chosen 
uniformly at random from Qn, converges (after rescaling) to a normal law as 
n — 00. 



1. Introduction 

Several enumeration problems on classes of labelled planar structures, e.g. la- 
belled planar graphs, were solved recently OHIHD]. These results were successfully 
used for ef&cient random generators for labelled planar graphs [T^ , based on Boltz- 
mann samplers [14] . In contrast, less is known about enumerative results on classes 
of unlabelled planar structures: the only classes treated so far are forests [25] and 
more recently outerplanar graphs [4]. 

In this paper we present a general framework to enumerate in a unified way a 
wide variety of labelled and unlabelled classes of graphs. Our main contribution is 
a universal method for rich classes of unlabelled graphs, the so-called "subcritical" 
classes of graphs, which is established through the novel singularity analysis of 
counting series and yields asymptotic estimates and limit laws for various graph 
parameters. In order to make this method more accessible and transparent, we 
include a brief analysis of the corresponding labelled classes, which was already 
carried out in [5T] . 

Another main contribution is asymptotic estimates and limit laws for various 
graph parameters of unlabelled series-parallel graphs. We study the class of unla- 
belled series-parallel graphs, firstly as an important subclass of planar graphs whose 
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asymptotic study has not been carried out so far and therefore it is interesting in 
its own right, and secondly as a concrete prototype-example to illustrate how our 
general method is applied. A graph is series-parallel (SP-graphs for conciseness) 
if its 2-connected components are obtained from a single edge by recursive sub- 
division of edges (series operation) and duplication of edges (parallel operation). 
Equivalently, SP-graphs can be defined in terms of minors as graphs which exclude 
K4 as a minor. Finally, a graph is series-parallel if and only if its tree-width is at 
most 2. Applying our general method, we show that the number of unlabelled 
SP-graphs on n vertices is asymptotically of the form 

g„ = cn-5/2p-" (1 + 0(1)), 

where p « 0.10655. Let G„ be a graph chosen uniformly at random among all 
unlabelled SP-graphs on n vertices. The random variable X„ counting the number 
of edges (blocks, or cut-vertices) in G„ features a central limit law 

%^&->A^(0,l), 
VVarX„ ^ ' 

where EX„ = fin + 0{l) and VarX„ = a'^n + 0{l) for computable constants fi and 
(7^. In addition, the random variable X!^ counting the number of vertices of degree 
k (for k fixed) in Gn satisfies a central limit law with mean EX^ = /ifeft + 0(1) 
and variance Yai — aln + 0{1) where fik and are computable constants. We 
also present a simple and general condition assuring > 0, which holds in a wide 
variety of graph families. 

Furthermore we show that the same subexponential term appears in other 

classes of graphs (in both labelled and unlabelled cases), which is, roughly speaking, 
inherited from a tree-like nature. This behaviour appears as a consequence of a 
subcritical composition scheme which appears in the specification of the counting 
series associated to connected graphs of the class. Such classes of graphs arising 
from subcritical composition scheme are called subcritical classes of graphs, whose 
formal definition is rather technical and therefore will be suspended to Sections 4 
and 5. 

We consider block-stable classes of graphs: we say a class G of graphs block-stable 
if and only if for each graph G G G, each of its 2-connected components (also called 
blocks) belongs to G- The class of cacti graphs, the class of outerplanar graphs, 
the class of SP-graphs, and other classes of graphs defined in terms of a class of 2- 
connected components are block-stable. Additionally, we consider classes of graphs 
which are defined by a finite set of 3-connected graphs. Observe that SP-graphs can 
be seen as graphs without 3-connected components. We show that the classes of 
cacti graphs, outerplanar graphs, and SP-graphs, are subcritical and prove that the 
asymptotic estimates and limit laws for graph parameters of subcritical classes of 
graphs follow the same asymptotic pattern and limit laws as the class of SP-graphs. 

The asymptotic study of subcritical classes of graphs consists of two steps: formal 
and analytic steps. The formal step consists in translating Tutte's seminal ideas on 
decomposing graphs into components of higher connectivity |71 1301 131j in terms of 
the decomposition grammar, which is comparable to the ones introduced in [5j (see 
also [I3)- This decomposition grammar translates combinatorial conditions into 
functional equations satisfied by the counting series of various classes of graphs. 
These counting series depend on the connectivity degree and the way how the 
graphs are rooted. In the analytic step, we extract singular expansions of the 



ASYMPTOTIC STUDY OF SUBCRITICAL GRAPH FAMILIES 



3 



counting series from the systems of functional equations. The main ingredient in 
this step is from [lOj . in which precise singular expansions are deduced for very 
general systems of functional equations. Finally, we derive asymptotic formulas 
from these singular expansions by extracting coefficients, based on the transfer 
theorems of singularity analysis [15l [16] . 

We also study natural parameters on a uniform random graph on n vertices cho- 
sen from a subcritical class, such as the number of edges, the number of blocks and 
the number of cut- vertices. These parameters satisfy normal limit laws, which come 
from the additivity behaviour of these parameters. We compute their expectation 
and variance, which characterise completely their limit distribution. These param- 
eters were studied in [21] for various labelled classes of graphs. In this paper we 
rediscover these results for labelled case and obtain new results for unlabelled case. 
Finally, applying the general techniques for systems of functional equations we de- 
duce the limit law for the degree distribution of a graph chosen uniformly at random 
among all the graphs on n vertices, by analysing in a unified way both labelled and 
unlabelled subcritical classes. The present work complement the work [3j [121 [E] 
for related problems on labelled graphs. Other non-additive parameters, such as 
the size of the largest block are not treated here. This parameter was studied in the 
labelled framework in [3T| and [55] . In the latter, the theory of Boltzmann samplers 
is applied. 

Outline of the paper. The paper is organised as follows. In Section [2] we 
introduce the notation and the terminology used in this paper. All the analytic 
machinery needed in order to deal with systems of functional equations, extraction 
of coefficients and limit laws is introduced in Section [3] In this section, we adapt 
the results from |10| to our context and recall the transfer theorems of singularity 
analysis from jl6) and other key results necessary to derive limit laws for graph 
parameters. We also introduce a novel result which assures the positivity of the 
variance under certain easy conditions. In Section [4] we obtain general results for 
the enumeration of labelled subcritical graph classes. These results are generalised 
to unlabelled subcritical graph classes in Section [sj These enumerative results we 
obtain are consequences of the general framework presented in Section [Sj Concrete 
examples are studied in Section [6] This section includes the analysis of unlabelled 
cacti graphs, unlabelled outerplanar graphs, and unlabelled SP-graphs. We also 
obtain general enumerative results for classes of graphs defined by a finite set of 
3-connected graphs. Limit laws are studied in Section [7] All parameters studied 
in this section give rise to normal distributed random variables, independently of 
the class (in either labelled or unlabelled setting). The degree distribution, a tech- 
nically more involved and interesting parameter, is studied in Section [8] Finally, 
the constant growth for unlabelled SP-graphs is computed in Section [9| using a 
numerical method. 

2. Graph classes, block-decomposition, and counting 

2.1. Combinatorial classes and counting series. As described in [31 [TH], a 

labelled combinatorial class is a set Q = U„>o Qn of objects such that Qn is finite for 
n > 0. Each object g in Qn has n labelled "atoms" (e.g., vertices of a graph) carrying 
distinct labels in {1, . . . ,n} and n is called the size of g. Two objects of G (of the 
same size n) are called isomorphic if one is defined from the other by relabelling. 
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It is always assumed that a combinatorial class is stable under relabelling, e.g., 
for a graph class Q we assume that a graph g is in Q if and only if all graphs 
isomorphic to g are also in Q. This way the symmetric group 6„ acts on Q„- for 
(T e (3„ and g Qn, o' ■ g has the same vertex and edge set as g, but each label 
i in g is replaced by a~^{i) in a ■ g and we write a ■ g = g. The set of objects of 
Gn considered up to isomorphism is denoted by Gn (in other words Gn = Gn/&n), 
and the combinatorial class G '■= U„ Gn is called the unlabelled combinatorial class 
associated with G- For counting purpose one classically considers the exponential 
generating function (shortly the EGF) in the labelled setting: 

^ — 'n>0 n'. 

and the ordinary generating function (OGF) in the unlabelled setting: 

^ — ^n>0 

For unlabelled enumeration, it proves convenient to consider a refinement of the 
OGF, called the cycle- index sum, a series in infinitely many variables si,S2,... 
defined as 

Ze(si,S2,...) := V -7 V 

('T,g)eS„xe„ 
"■9=9 

where :— s^^Sj^ • • • sj^" is the weight- monomial of a permutation a of cycle type 
l'=i2'^^ • • • n'^" (that is, Ci cycles of length i for z = 1, . . . , n). The OGF is recovered 
from a specialisation of the cycle index sum |22] by replacing Si by for each i > 1: 

G{z)^Zg{z,z\z\...). 

We will consider classes of graphs of various types depending on whether one marks 
vertices or not. All graphs are assumed to be simple (no loops nor multiple edges) 
and are labelled at vertices. A (vertex-) rooted graph is a graph with a distinguished 
(labelled) vertex. A derived graph or pointed graph is a graph where one vertex 
is distinguished but not labelled (the other n — 1 vertices have distinct labels in 
1}). Isomorphisms between two pointed graphs (or between two derived 
graphs) have to respect the distinguished vertex. 

Given a graph class G, the rooted class G' is the class of rooted graphs from G, 
and the derived class G' is the class of derived graphs from G] since — n\Gn\ 

and \Gn\ n\Gn\, we have respectively G'{z) = j^Giz) and G*{z) = zG'{z). 

2.2. Block-decomposition of a graph. For fc > 0, a graph is fc-connected if one 
needs to delete at least k vertices to disconnect it. Obviously, a graph G is a set 
of its connected components. For the decomposition from connected graphs into 
2-connected graphs we use the block structure of a connected graph. A block of a 
graph G is a maximal 2-connected induced subgraph of G. We say a vertex of G 
is incident to a block B of G if it belongs to B. The block structure of G yields a 
bipartite tree with the vertex set consisting of two types of nodes, i.e. cut- vertices 
and blocks of G, and the edge set describing the incidences between the cut- vertices 
and blocks of G. This suggests a natural decomposition of connected graphs into 
2-connected graphs and this holds also for rooted graphs. The root-vertex of a 
rooted graph G is incident to a set of blocks and to each non-root vertex on these 
blocks is attached a rooted connected graph. In other words, a rooted connected 
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graph rooted at ! w is uniquely obtained as follows: take a set of derived 2-connected 
graphs and merge them at their pointed (distinguished but not labelled) vertices so 
that V is incident to these derived 2-connected graphs, then replace each non-root 
vertex w in these blocks by a rooted connected graph rooted at w (which is allowed 
to consist of a single vertex and in this case it has no effect). 

Through the entire paper, given a class of graphs Q, we denote by C (resp. B) the 
subfamily of connected (resp. 2-connected) graphs in Q. In the language of symbolic 
combinatorics from ^ [TB] , the block-decomposition described above translates into 
the fundamental equations: 

(1) g = Set(C), 

(2) C* =Z-Set(S'oC*), 

where the factor Z in the last equation takes account of the root vertex (which 
is distinguished and labelled), the symbol • denotes the partitional product on 
combinatorial classes, and the symbol o denotes substitution at an atom (see [2] 
for definitions) . As shown in Table [l] there is a well-known dictionary ^ [T^ , both 
in the labelled and in the unlabelled setting, that translates equations relating 
combinatorial classes into equations relating the associated counting series. 



Construction 


Class 


Labelled setting 


Unlabelled setting 


Sum 


C 


= A + B 


C{z) = Aiz) + Biz) 


C{z) = A{z) + B{z) 


Product 


C 


^A-B 


C{z) ^ A{z) ■ B{z) 


C{z) ^ A{z) ■ B{z) 


Set 


c 


= Sct(6) 


C{z) = cxp{B{z)) 


C(z)=exp(E,>a^(^0) 


Substitution 


c 


= AoB 


C{z)=A{B{z)) 


C(z) =Z^(6(z),^(z2),...) 



Table 1. The dictionary that translates combinatorial construc- 



tions into operations on counting series. 



A graph class Q is called block-stable if it contains the link-graph £, which is a 
graph with one edge together with its two (labelled) end vertices, and satisfies the 
property that a graph G belongs to Q if and only if all the blocks of G belong to Q. 
Block-stable classes include classes of graph specified by a finite list of forbidden 
minors that are all 2-connected, for instance, planar graphs {Forhid{K5, Ks^^)), 
series-parallel graphs (Forbid(if4)), and outerplanar graphs (Forbid(fir4, if3,2))- For 
a block-stable graph class, ([T]) and ^ translates into equations of EGFs in the 
labelled setting: 

Giz) ^ exp(C(z)), 
C'{z) = zexp(S'(C*(z))), 
and into equations of OGFs in the unlabelled setting: 

g{z) = exp icV)) , 

C^(z) = zexp (j2,^^lZB'iC'{z^),C^'{z''),C^'iz'^),...)). 
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A refined version of tlic last equation will turn out to be useful later, which expresses 
([2]) in terms of the cycle index sum: 

(3) Zc'{si,S2, . . .) = Si exp (X!j>i l^B'{Zc-{Si, S2i, ■ ■ .), ^C'(s2i, S4i, ■■■),■■ ■)) • 

With these systems of equations of EGFs, OGFs and the cycle index sum, 
we will perform singularity analysis under certain general conditions - so-called 
"subcriticallity" conditions - in order to get asymptotic results (asymptotic enu- 
meration and limit laws for graph parameters). The tools for singularity analysis 
are described in the next section. 

3. Tools for the asymptotic analysis 

The purpose of this section is to collect several facts on solutions of functional 
equations. Most of the following properties (and proofs) can be found in pT[ ITB] . 
which serve as general references for this subject. For our purposes we, however, 
need to adjust several points with additional properties. 

3.1. Singular expansions of multivariate series. We consider a power series 
of the form y = y{z, ui, . . . , Vk) in general with nonnegative coefficients, where z is 
singled out as the primary variable and vi, . . . ,Vk are secondary variables (possibly 
there is none). From now on we use the abbreviation v for (wi, . . . ,Vk). A valuation 
Vq of V is called admissible if all components of Vq are positive and if y(z, Vq) is a 
valid power series in z, i.e., [z"]j/(z, Vq) < 00 for each n > 0. 

Consider a fixed positive valuation (zqjVo) of (z,v). Then y is said to have a 
square-root expansion around (zo,Vo) if vg is an admissible valuation, zq is the 
radius of convergence of 2/(2:, vq), and the expansion 

y{z,v) = a(z,v) - b{z,v)^/l - z/p(v) 

holds in a neighbourhood of (zo,Vo) (except in the part where 1 — z/p{\') G M^), 
where the functions a(z,v) and &(z,v) are analytic at (zo,vo), fo(2;o,vo) > 0, p(v) 
is analytic at vq, and p(vo) = zq. The function /o(v) in the expansion is called the 
singularity function of y relative to z. 

Moreover, we will be particularly interested in functions j/(z,v) with such a 
singular behaviour, where /9(v) is the only singularity on the circle \z\ — \p{'v)\ (if v 
varies in a suitable neighbourhood of vo) and y[z, v) can be analytically continued 
to the region {2; G C : |z| < |p(v)| -I- e, 1 - zj p{w) ^ M^} (for some e > that 
is uniform in this neighbourhood of vq). In this case one can use an asymptotic 



transfer principle by Flajolet and Odlyzko jTS] (see also Section 3.31 to obtain 
asymptotics for the coefficient [z"] i/(z, v) of the form 

y{z, V) = ^^^^ n^l^ p(v)-" (H- O (n-^)) . 

Similarly, y is said to have a singular expansion of order 3/2 at (zo,vo) if the 
expansion is of the form 

y(z, v) = a(z, v) + &(z, v) • (1 - z/p(v))3/2, 

with otherwise the same conditions as for square-root expansions. 

It is not difficult to show that if a function ?/(z, v) admits a square-root expansion 
around (zq, Vq), then the function ?/(x, v) dx admits a singular expansion of order 
3/2 at (zo,vo) (see [H]). 
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3.2. Singularity analysis of systems of functional equations. Consider a 
vector y = {yi, . . . ,yj.) of formal power series in the formal variables z,v, with 
V = (vi, . . . , Vk) which is a solution of an equation-system of the form 

(E): y = F(y;z,v), 

where F(y; z, v) = (i^i(y; z, v), . . . , Fr{y; z, v)) are power series with nonnegative 
coefficients. In order to avoid trivial situations as in equations like y ~ zy"^, where 
y{z) = is the only solution, we will assume that y = is not a solution. We call 
such a system positive system. 

The Jacobian matrix Jacp of the system is the rxr matrix whose (i, j)-coefiicient 
is dFj/dyi. The singularity equation associated with (E) is 

(S) : = Det (Ir - Jacp) , 

where Ir is the rxr identity matrix. The singularity system of (E) is the system 
{(E), (S)}, which has r-t- 1 equations (the first k ones for (E), the last one for (S)). 

The dependency graph of F is the directed graph on V = {1, . . . ,r} such that 
there is an edge from i to j if and only if the jth component Fj of F really involves 
yi, that is, the power series dFj/dyi is not 0. The system (E) is called strongly 
recursive if the dependency graph of F is strongly connected (which means that 
every pair of vertices (wi,W2) is linked by a directed path). 

Informally this condition says that no subsystem of (E) can be solved prior to 
the whole system of equations. An equivalent condition is that the correspond- 
ing adjacency matrix and, thus, the Jacobian matrix Jacp is irreducible. The 
most important property of irreducible matrices A with non-negative entries is 
the Perron- Frobenius theorem (see [24] ) saying that there is a unique positive and 
simple eigenvalue Xmax = Amaa;(A) with the property that all other eigenvalues A 
satisfy |A| < Xmax- This unique positive eigenvalue is a strictly increasing func- 
tion of the entries of the non-negative matrix. More precisely, if A = (a.^) and 
A' = (a'y ) are different irreducible non-negative matrices with aij < a'^^ (for all 
i,j) then X^axi-^) < -^maa;(A'). Moreovcr, every principal submatrix has a smaller 
dominant eigenvalue. 

Consider an admissible valuation vq of v and assume that the radius of conver- 
gence ^0 of z — >■ yi{z, Vq) is finite and strictly positive. If the system (E) is strongly 
recursive then it is easily shown that, for each i g {1, . . . zq is the radius of 
convergence of yi (due to the dependency graph being strongly connected) and 
yi(z,Vo) converges to a finite constant > as z ^ zq (due to the fact that F is 
nonlinear according to y). 

Note also that (S) says that A = 1 is an eigenvalue of Jacp. We observe that 
if y = y(-2,vo) is an analytic solution of a strongly recursive system (E) that is 
singular at z = zq and yo = y(^o,vo) is a finite vector, then (S) is satisfied for 
(2/;z,v) — (?/o;2o,vo) provided that (yo;zo,vo) is an inner point of the region of 
convergence of F. (This explains the term singularity system.) However, in order to 
obtain the radius of convergence (in positive and analytically well-founded systems) 
we will need the condition that A„iaK(JacF) = 1 (see [T]). 

Definition 1. A system (E) is called analytically well-founded at a fixed positive 
valuation Vq if the following conditions are satisfied: 

(1) The valuation Vq is admissible for the system of power series F(y;z,Vo); 
we have [y^z"] i^i(y;z,vo) > for all i, and F(0;z,vo) 7^ 0. 
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(2) The function F is not affine in y and depends on z, that is, there are i and 
j in {1, . . . ,r} such that d'^'F / dyidy.j ^ and d¥ jdz ^ 0, 

(3) There exist > and a positive vector yg for which (yo; zq, Vq) is an inner 
point of the region of convergence of Fi{y; z, v) (for i S {1, . . . , r}j and (E) 
and (S) are satisfied for (y;z,v) = (yo;^OiVo) and that all eigenvalues A 
of the Jacobian matrix 3a.Cp{yQ; zq,vo) satisfy |A| < 1. 

It is not clear that the singularity system (E), (S) has a proper solution. However, 
if (E) is a positive and strongly recursive system and if F(y; z, v) is a vector of entire 
functions in y and z (for v in a neighbourhood of vq) then it is always solvable. 
In particular it is enough to consider the singularity of the power series solution of 
(E). We also get the property Xmaxi^acpiyo; zo,vo j) = 1. 

Note that the radius of convergence zq of the solution y of a positive and strongly 
recursive system (E) is always finite. Furthermore, y{zo) is finite, too, if F is not 
affine in y. 

The following theorem contains the main properties of systems of equations (E) 
that we will use in the sequel (for a proof see [llj). 

Theorem 2. Suppose that an equation (E) is strongly recursive and analytically 
well-founded at a valuation Vq. Then there is a unique vector of power series 
y = y(z, v) in the variables z,v that satisfies (E). Furthermore, the components of 
y have non-negative coefficients [z"] yi(z, Vq) (for i € {1, . . . ,r}) and a square-root 
expansions around (zo,Vo). 

Moreover, if [z"] yi{z, Vq) > for all n > Hq then Zg is the only singularity on the 
radius of convergence |z| — zq and all components y can be analytically continued 
to the region D = {z & C : \z\ < |/9(v)| + e, 1 — z/p{v) ^ where e > is 

uniform for v in some neighbourhood of vg . 

The condition [z"] 2/i(z, vg) > (for n > rtg) is usually verified by using a 
combinatorial interpretation of the coefficients. In the case of a single equation 
y = F{y;z,v) = J2n m k ^n,m,kz"'y"^v'' it is also possible to check this with the 
help of conditions on the coefficients an^m,k- For example, if there exist m > 1, 
ni,n2,n3 and fci,fc2,^3 with am,m,ki 7^ 0, On^^^M 0, an3,m,k3 7^ such that 
n2 — ni and — ni are coprime or if there are ni,n2, mi > 1, > Ij ^11^2 with 
a,u,mi-i,fci 7^ 0, an^^jn2-iM such that ni{m2 - 1) -n2(mi - 1) = 1 then it also 
follows that [z"] y{z,vo) > (for n > Tig) - compare it with and the methods 
used in the proof of Lemma [4j 

3.3. Transfer theorems of singularity analysis and central limit theorems. 

As detailed in the book by Flajolct and Sedgewick the singular behaviour of 
a counting series can be translated to an asymptotic estimate of the counting se- 
quence, by coefficient extraction (which is done via contour integrals in the complex 
plane). 

Let Zg be a non-zero complex number, and e and S positive (real) numbers. Then 
the region 

A = A(zg,e,(5) = {z £ C : \z\ < zo + e, | arg(z/zg - 1)| > 5} 

is called a A-region. The basic observation (see [15]) is that a singular expansion 
around the singularity Zg that is uniform in a A-region transfers directly to an 
asymptotic expansion for the coefficients. Suppose that a function y{z) is analytic 
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in a region A-region A(zo, e, S) and satisfies 

y{z) = C (1 - z/zof + O ((1 - z/zof) , z G A(zo, e, S), 
where /? > a and a is a non-negative integer. Then we have 

(4) [z"] y{z) = Cj^zo- + O ^-^-i}) . 

It is an important additional observation that the imphcit constants are also effec- 
tive which means that the O-constant in the expansion of y{x) provides explicitly 
an O-constant for the expansion for [z"] y{z), and that the same statement is true 
if we change the O-constants by o-constants. See [TS] for details. In particular it 
follows that singular expansions that are uniform in some parameter also translate 
into asymptotic expansions of the form Q with a uniform error term. In particular 
it applies for functions y(z,v) with square-root expansion around (zq^vq) or with 
singular expansion of order 3/2, provided that they can be analytically continued 
to a A-region. 

Next we restrict ourselves to univariate w € C and are interested in bivariate 
asymptotic expansions of the coefficients [z"v'^]y{z,v) when y has a square-root 
expansion. We introduce the function fi defined as 

vp'{v) 
p(v) 

We call it regular in a closed interval V of the positive real line if is strictly in- 
creasing on V . In this case v = p,~^ denotes the inverse function of /i. Furthermore, 
we set (t{v) — ^Jvpljv) which is positive in the regular case. If we additionally as- 
sume that zq = p(wo) is the only singularity of y(z, v) for |z| — zq and \v\ — vq then 
we have uniformly for m/n € i^(V) 



p{v{m/n))-''v{m/n)-"' (l + O {n-^)) 



(5) [z"z;™]y(z,«) 

b{p{v{m/n))^ v{m/n)) 
2-\/2 7r {n/m)a{u{m/n)) mn 

(For a proof we refer to ?]). Again it is clear that this has a direct analogue for 
functions with a singular expansion of order 3/2. 

It is relatively easy to check the condition that zq = p{vo) is the only singularity 
of y{z, v) for |z| = zq and \v\ — vq and also that p is regular - compare it with the 
remark following the proof of Lemma |4] 

Suppose that m = pn + 0{\/n). The asymptotic expansion ([5| for the coefficient 
[z"w™] y(z, v) behaves locally like 

(6) [z%™] 2/(z, v) = J^^^WlA) e-(™-A^(i)»)V(2-(i)^n)^(i)-« ( 1 + o(„-i/2)) 

2v7rcr(l)7i V / 

which suggests that there is a central limit theorem behind. Actually this is true. 

Let Xn be a random variable with probability distribution P(A'„ = m) = 
[z"i;™] y{z, ?j)/[z"] j/(z, 1), then the asymptotic expansion ^ is indeed a local limit 
theorem for X„. In general, there is a combinatorial central limit theorem. For the 
sake of brevity we do not list all possible versions but only for a single equation 
and we comment on systems of equations. 
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Theorem 3. Suppose that Xn is a sequence of random variables whose probability 
generating function is given by 

" = 7 1 — , 

[x"] yiz, 1) 

where y{z, v) is a power series that is the (analytic) solution of the functional equa- 
tion y = F(y;z,v), where F{y;z,v) satisfies the assumptions of Theorem^ In 
particular, let z = zq > and y = yo > be the proper solution of the system of 
equations y — F{y; z, 1), 1 — Fy{y; z, 1) Set 

F, 
zqF^ 

1 

ZoFfFyy 

+ {PyyPzz ~ Pyz)^ i 

where all partial derivatives are evaluated at the point (yo; zq,1). Then the asymp- 
totic mean and variance of X„ satisfy 

EX„ =/in + 0(l) and VarX„ = cr^^ + 0(1) 

and if > 



^ l^ ^ l^ ^" _ ( {-^yy-^vv ^yv) '^-^z^vi-^yy-^zv ^yx^yv) 



Xn — E Xn 

y/YarXn 



N{0,1). 



Note that u, = is the same as the other n defined above as u(v) = — ""^ 
when v = \. Similarly we have tr^ = cr^(l). 

The situation for a system of equations is quite similar (even if we additionally 
consider a random vector X„ instead of a random variable Xn)- Suppose that y = 
(j/i, . . . ^yr) is the solution of a system of equations (E) and that the assumptions 
of Theorem [2] are satisfied. Furthermore set y(z,v) = iJ(y(z, v); z, v) for a power 
series F with non-negative coefficients, for which (yo;2;ojl) is inner point of the 
region of convergence and we have i/y(yo;zo,l) ^ 0. Then the random vector 
X„ = (Xi;„, . . . , Xfc;„) with probability generating function 

„ Y „ X, X, \z^^] viz, v) 

1 ^ [z"] y(z, 1) 

is asymptotically normal with asymptotic mean E X„ = n/x + 0(1) and covariance 
matrix CovX„ — nS + 0(1), where 

1 b'^Fv(yo;zo,l) 
^ zobTF,(yo;zo,l)' 

in which b is (up to scaling) the unique positive left eigenvector of Jacp, and S 
is a positive semi-definite matrix which can be computed with the help of second 
derivatives (for details see [H]). In many applications b appears to be (1, . . . , 1)'^, 
which is due to the special structure of the systems of equations. The source 
of the central limit theorem is actually a singular expansion with singular term 



1 ^ with a and thus a central limit theorem with the same mean and 

^For convenience we use the notation Fy to denote the partial derivative dF/dy. 
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variance follows also for generating functions given by / y{z,^) — J F{y{z,v); z,v), 
which have the same singularity, but of order | . 

Finally we comment on the positivity of cr^ = cr^( 1) in the case of a single func- 
tional equation y = F{y;z,v). (Equivalently this concerns the question whether 
^(v) = vF^{y(p{v),v); p{v),v)/{p{v)F;,{y{p{v),v); p{v),v)) is regular in a neigh- 
bourhood of u = 1.) 

Lemma 4. Let y — F{y,z,v) = X^n m fe ^".'".'=•^"2/'"''^'^ '^'^ analytically well 
founded equation for the valuation Vq — 1 as given in Theorem Suppose that 
there are three triples {nj,mj,kj) , j = 1,2,3, of integers with ruj > 0, j — 1,2,3, 
and 

ni mi — 1 fci 

712 "T-2 - 1 k2 ^0 

713 - 1 fcs 

such that an-^mj,kj 7^ 0; i = 1- 2,3. Then > 0. 

Proof. Let x — p{v) be the solution of the singular system y — F{y;z,v), 1 = 
Fy{y; z,v). We will first show that |/5(e**)| > p{l) for real t ^ that are sufficiently 
small. This property will be then used to prove that cr^ > 0. 

First it is clear that |/5(e**)| > p{l) for all real t for which p(e**) exists. For, if 
|p(e'*)| < p(l) then we would have |2/(p(e**, e**)| < j/(p(l), 1) and also 

(8) \Fyiyipie^\ e^*); p(e^*), e»*)| < Fyiyip{l), 1); p(l), 1) = 1 

which is a contradiction. (Note that we have used here the assumptions Fyy ^ 
and F„ ^ 0.) 

Now assume that |p(e**)| < p(l) for some real number t. Then an inequality 
similar to ([s]) implies that for all n, to, k 

^o,n,m,k^ y V = rna^ jjij^ZQl/Q , 

where we used the abbreviations zq = p{l), yo = 2/(p(l),l), z = p(e'*), y = 
y(/9(e**, e**), and v = e'*. In particular, it follows that z"^^ y'^^^^v^^ = z^^y^^ ^ for 
j — 1, 2, 3. Hence, if we set z = zoe"", y — y^e'^'^ (and v — e'*) it follows that 

Ujr + {mj — \)s + kjt ~ 2Trlj (j = 1, 2, 3) 

for some integers Ij, j — 1, 2, 3. This is a regular system and implies that there is a 
(unique) solution of the form r = 27rLi/M, s = 2ttL2/M, t = (for integers 

Li,L2,L3,M). Hence, if i ^ is sufficiently close to then |/c(e**)| > p{l). 
Next consider the Taylor series of the function 

oo 

g{t) ^ log p{e*)=Y,'^t' ■ 

3=0 ^' 

By definition we have kq — logp(l), ki — p,, and K2 = . Note that this represen- 
tation and the general property |p(e'*)| > p(l) implies that k.2 > 0. Suppose that 
K2 = and let i'o > 3 be the smallest integer with k^^ ^ 0. 

We use now the fact that the assumptions of Theorem |3] imply that 

Ee-" = (g*^)"" (1 + 0(0). 
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(This follows from the singular expansion of the solution y(z, v) and the asymptotic 
transfer results - compare it with Q from above). Hence, by using the Taylor 
expansion of g{t) it follows that 

This means that the sequence of random variables Yn — (X„ — fj,n)n~^^^° converges 
weakly (and we have convergence of all moments) to a random variable Y with 
Laplace transform E e*^ = e~'^*o* °/^o\ However, such a random variable that has 
a non-zero £o-th moment but zero variance does not exist. Hence, we finally have 
proved > 0. □ 

A slight variation of the above proof shows that if there are three triples {uj , rrij — 
l,kj), j = 1,2,3, with determinant ±1 then |p(e'*)| > p(l) for aU t ^ 2nZ. This 
shows that zq = p(l) is the only singularity of y(z,v) for \z\ — zq and \v\ = 1. This 
assumption can be used to obtain bivariate asymptotic of the form ([5|. 

4. SUBCRITICAL GRAPH CLASSES: THE LABELLED CASE 

In this section Q denotes always a block-stable class of labelled graphs and C 
(resp. B) its subclass consisting of connected (resp. 2-connected) graphs. 

4.1. Definition of subcriticality. Recall from Section [2^2] that the EGFs of the 
block-stable class satisfy 

g{z) = cxp(C(z)), C'iz) = zcxp(S'(C*(z))). 

Given a series g{y) it is easy to show that there is a unique series f{z) that is a 
solution of the equation 

(9) /(z) = zexp(5(/(z))). 

In addition f{z) has nonnegative coefficients if g{y) has nonnegative coefficients. 
Note that for a block-stable graph class, the solution of Q is /(z) = C*{z) when 
g{y) is taken as B'{y). 

Definition 5. Let g{y) he a series with non-negative coefficients such that g(0) = 0. 
Let f{z) be the unique solution of Let p and rj be the radii of convergence of 
z I—)- /(z) and y i-> g{y)- Then the pair (f{z),g(y)) is called subcritical if f{p) < rj. 

A block- stable graph class Q with C the connected subclass and B the 2-connected 
subclass is called subcritical if the pair {C*{z),B' [y)) is subcritical. 

Note that the singularity system for ([9| is 

y = zexp(g(y)), y g' {y) = 1. 

In particular, for the pair {C [z) , B' {y)) the latter equation rewrites to yB"{y) = 1. 
Hence, we have subcritiallity if and only if r]B"{r]) > 1, compare with 

In what follows we will also consider functions /(z, u), g{y, v) with an additional 
parameter v. For example, suppose that we are dealing with a bivariate generating 
function where the exponent of v counts the number of edges and the exponent 
of z (or y) counts the number of vertices. Suppose further that we already know 
that the pair {f{z,l),g{y,l)) is subcritical. What can we say then for the pair 
{f{z,v),g{y,v)) if v is sufficiently close to 1? Is there some stability of the subcriti- 
callity? Actually there is if the the parameter that is counted by the exponent of v 
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is at has a linear worst case behaviour in the exponent n of z (or y). The essential 
consequence of the following lemma is that radius of convergence z i— > f{z,v) and 
y 5(2/1 '^): respectively, is continuous at u = 1. Hence we have f{p{v),v) < rj{v) 
if V is real and sufficiently close to 

Lemma 6. Let A{z, v) — ^„ j,>q o„^fcz"u'^ be a power series with non-negative 
coefficients On^k with the property that an,k — for k > Cn for some constant 
C > 0. Let R{v) denote the radius of convergence of the mapping z 1— > A{z^v). 
Then we have for real v > 

iT!(l)min{l,w"<^} < R{v) < i?(l) max{l, u*^}. 
Proof. If v > 1 then 



an,k < ^ an^kv'' < ^ an,k v'^" 

k k \ k J 



and consequently R{v) > R{l)v . Similarly we argue for < v < 1. □ 

4.2. Asymptotic estimate for a subcritical graph classes. We start with 
a quick analysis of subcritical graph classes and derive their asymptotic number, 
compare it with [3j. 

Lemma 7. Let Q be a labelled subcritical block-stable graph class with C the con- 
nected subclass and B the 2-connected subclass. Then C*{z) has a square-root sin- 
gular expansion around its radius of convergence p. Furthermore, if [z^^]C*{z) > 
for n > no then p is the only singularity on the circle \z\ — r and C*{z) can be 
continued analytically to the region D = {z S C : |z| < p + e, 1 — z/ p ^ M~} for 
some e > 0. 

Proof. The function y = C*(z) is a solution of 

y = F{y; z), with F[y, z) = z exp(S'(j/)). 

Let p and 77 be respectively the radii of convergence of C'{z) and of B'{y), and let 
T := C'{p). Since G is subcritical, we have t < rj, hence F{y; z) is analytic at (t, p). 
We conclude from Theorem [2] that C*{z) has a square-root expansion at p and that 
C*{z) can be continued analytically to D. □ 

Theorem 8. Let Q be a subcritical block-stable graph class with the property that 
[z'"']C*{z) > for n > uq. Then there exist constants 7 > e « 2.71828 and c > 
.such that 

(10) [z'']g{z) = cn-'^^^Y\l + o{l)) as 71 ^00. 

Proof. The function C{z) satisfies 

C(z)=^'r(i)y, 

hence C{z) has a singular expansion of order 3/2 at p. Since Giz) = exp(C(z)) and 
exp is analytic everywhere (in particular at C{p)), we conclude that Giz) has also a 



^Note that it is sufficient to consider real v if we are just interested into (global) central 
limit theorem and asymptotics for moments. Namely, in order to prove a theorem of the type 
of Theorem [s] one can work with the help of the Laplace transform E e'^" that is encoded by 
]ggtx„ ^ [z"]A(z,e'^)/lz"]A{z, 1) when F{X„ = k) = a„,k/a.n. 
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singular expansion of order 3/2 at p. The transfer theorems of singularity analysis 
(see Section 3.3 ) yield an estimate of the form ( 10 1, where 7 = 1/ p. See also [3]. □ 

4.3. Sufficient condition for subcriticality. The following lemma gives a simple 
sufficient condition for subcriticality: 

Lemma 9. Let g{y) he a series with non-negative coefficients and positive radius 
of convergence rj such that g' {y) — > 00 as y ^ r]~ . Let f{z) be the unique solution 
of ([9|. Then the pair {f{z),g{y)) is subcritical. 

Proof Let p be the radius of convergence of f{z) and t := /(/?). From the definition 
of the subcriticality we need to show t < rj. 

Assume t > rj. Then, by continuity of /(z), there exists < zq < p such that 
/(zq) — rj. Since /(z) is regular at zq with positive derivative and since g{y) is 
singular at 77, the function g{f{z)) must be singular at zq. Hence z exjp(g{f{z))) 
must also be singular at zq, in contradiction to the fact that /(z) is regular at zq. 
Hence t < r/. 

Assume now that t = 77. Differentiating ([9|, we obtain 

nz) = ^+g'{fiz))riz)f{z), 

which imphes /'(z) > g' {f{z))f' {z)f{z), for z £ (0,p). Taking y = /(z), it sim- 
plifies to g'{y) < 1/y for y E (0,7]). This contradicts the fact that g'{y) — cx) as 
y ^ rj. □ 

Note that the pair {f{z),g{y)) is subcritical if g{y) has a square-root singular 
expansion at rj. Therefore, a block-stable graph class G is subcritical if the EGF of 
B' admits a square-root singular expansion. 

5. Subcritical classes: the unlabelled case 

In this section Q denotes a block-stable class of unlabelled graphs and C (resp. 
B) its subclass consisting of connected (resp. 2-connected) unlabelled graphs in Q. 

5.1. Definition of subcriticality in the unlabelled case. We have seen in 
Section 12.21 that a block-stable class satisfies 



g{z) = cxp (^^^^ iC(zO 



C'iz) = z exp \ZB'iC'{z^),C'iz'^),C'{z'l, ■ ■ 

The second equation can be rewritten as follows: 

(11) /(z) = zexp(5(/(z),z)-fA(z)), 
where 

(12) /(z) := C'iz) 

(13) giy,z) := ZB'{yJiz^)Jiz'),...) 

(14) Aiz) := ^^^^izB'(/(^^),/(^''X---)- 

Note that, given a bivariate series g{y,z) and a univariate series A{z), there is 
a unique series /(z) that is a solution of (11) (because the coefficients of /(z) are 
determined uniquely iteratively) and /(z) has nonnegative coefficients if g{y, z) and 
A{z) have nonnegative coefficients. 
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Definition 10. Let g{y, z), A{z) be series with nonnegative coefficients, and let 



f{z) be the unique solution of (111 and p the radius of convergence of f{z). Then 
the triple {f{z),g{y,z),A{z)) is called subcritical if 

(i) p is non-zero, 

(ii) g{y,z) is analytic at {f{p),p), and 

(iii) the radius of convergence of A(z) is larger than p. 

An unlabelled block-stable graph class Q is called subcritical if 



(a) the triple {f{z),g[y,z),A{z)) with f{z), g{y,z), A{z) defined as in (12) 



( 14 ) is subcritical, and 
(b) the radius of convergence of the series q{z) := Zc{0, z'^ , z^ , . . .) is strictly 
larger than p. 

Note that for any block-stable graph class G, the class C* of rooted connected 
graphs from Q dominates coefficient-wise the class of unlabelled rooted non-plane 
trees, whose coefficients grow exponentially; hence, p < p* ~ 0.33832 (with p* the 
radius of convergence of unlabelled forests). Furthermore, we also have stability 
of subcriticallity when we vary an additional variable v locally around 1 if the 
parameter that is counted by the exponent of v has at most linear worst case 
behaviour. 

5.2. Asymptotic estimate for a subcritical class. 

Lemma 11. Let Q be an unlabelled subcritical block-stable graph class with C the 
connected subclass and B the 2-connected subclass. Let p be the radius of conver- 
gence of f{z) := C*{z). Then f{z) has a square-root singular expansion around p, 
and {y, z) — {f{p), p) is a solution of the singular system 

y = z exp{g{y,z) -\- A{z)), l^ygy{y,z), 



with g{y,z) and A{z) defined in (13) and (14). Furthermore, if [z"]C*(z) > for 
n > uq then p is the only singularity on the circe \z\ — r and we C*{z) can be 
continued analytically to the region D = {z G C : \z\ < p -\- e, 1 — z/p ^ K~} for 
some e > 0. 

Proof. Recall that the function f{z) ~ C*{z) is a solution of 

y = F{y;z) ^ zexp{g{y,z) + A{z)). 

Note that h{z) = Zigr{f{z),f{z'^),...) is bounded coefficient-wise above by f{z) 
and hence the singularity of h{z) is larger than p. Since p G (0, 1) (by the remark 
just after Definition 10), the function A{z) = X]i>2^(^*)/* analytic at p. By 



definition of subcriticality also, the function g{y,z) is analytic at {f{p),p). Hence 
F{y] z) is analytic at {f{p), p). Since the system is clearly strongly recursive and the 
function f{z) aperiodic, we conclude from Theorem [2] that f[z) has a square-root 
expansion at p. □ 

Lemma 12. Let Q be an unlabelled subcritical block-stable graph class with C the 
connected subclass and B the 2-connected subclass. Let p be the radius of conver- 
gence of f{z) :— C*{z). Define R{s,z) :— Zc'{s, z'^ , z^ , . . .). Then R{s,z) has a 
square-root singular expansion around {p,p), and the singularity function ^(z) of 
s I— >■ R{s,z) has a negative derivative at p. 
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Proof. The bivariate series R{s,z) is a refinement of C*(z), since C*{z) = R{z,z). 
Tfie equation ([S]) implies tliat y = F{y;z,s) := s exp{g{y, z) + A{z)), with g{y,z) 



and ^(z) defined in (13) and ([14|. The singular system for R{s,z) is 



y = s exp(.g(y, z) +^(z)), 1 = ygy{y,z). 



This is the same as the singular system of f{z) (given in Lemma 11 1 except that the 
variable z on the left-hand side of exp is replaced by the variable s. By Lemma [Tl] 
(y, z) — {f(p),p) is a solution of the singular system of f{z), hence clearly (y; z, s) = 
(/(p); Pi P) is a solution of the singular system of i?(s, z), and z, s) is analytic 
at {f{p);p,p), since g{y,z) is analytic at {f{p),p). Thus, Theorem [2] ensures that 
R{s,z) has a square-root singular expansion at {p,p). In addition, the singularity 
function ^(z) has a negative derivative, since F(y; z, s) depends only on z. □ 

Theorem 13. Let Q he an unlabelled subcritical block-stable graph class such that 
[z"]C*(z) > for n> riQ. Then there exist constants c > and 7 such that 

[z"']g{z)=cn-^/^Y'{l + o{l)) asn^oo 

for 7 > 7* ~ 2.95576, where 7* is the growth rate of unlabelled forests. 

Proof. First we show that C(z) has a singular expansion of order 3/2 at p. Define 
Q{s,z) :— Zc{s, z"^ , z^ , . . .) (note that C(z) = Q{z,z)). The general relation Z^^i = 
ij^Zj^ ensures that R(s,z) = sQs{s,z), hence 

Q{s,z) = Q{0,z)+ / R{w,z) 



JO w 

,2 ,3 



The term Q(0, z) = Zc{0, z , z , . ■ .) — q{z) is analytic at p, by definition of sub- 
criticality. Since R{s, z) has a square-root expansion at (p, p), the integral term has 



a singular expansion of order 3/2 at (p, p) (see Section 3.1) of the form 

Q(s, z) = a(s, z) &(s, -2) • (1 - sM-z))^/'. 
Therefore, C(z) = Q{z, z) has a singular expansion of the form 

C(z) = a(z, z) + biz, z) ■ ((p(z) - z)lp{z)f\ 

Since ^(z) has a negative derivative at p and ^(p) = p, there exists a function A(z) 
analytic and nonzero at p such that p(z) — z = (1 — ^^/p) ■ A(z). We conclude that 
C(z) has a singular expansion of order 3/2, of the form 

C{z)^a{z)+P{z)-{l^zlpf\ 

with a{z) = a{z, z) and (3{z) — b(z, z) ■ (A(z)/p(z))'^^^. 
Recall that Q{z) and C(z) are related by 

g{z) = exp(C(z) + S(z)), with £;(z) 7^(^')- 

Since i:^(z) is analytic at p, the singular expansion of order 3/2 at p for C(z) yields 
also a singular expansion of order 3/2 at p for Q{z). The transfer theorems of 
singularity analysis then yield the estimate for [z"]CJ(z). □ 
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5.3. Sufficient conditions for subcriticality. Similarly as in the labelled case, 
we provide a list of conditions that implies subcriticality, but will be convenient to 
check on examples (see Section |6] for the application) : 

Lemma 14. Let Q he an unlabelled block-stable graph class with C and B the con- 



nected and 2-connected subclasses. Let f(z),g{y,z) be defined as (12) and (13) 
and let p be the radius of convergence of f{z). For z > let rj{z) be the radius of 
convergence of y g{y,z). Assume that 

(1) there exist constants c and 7 > such that [z'^]C* < C7", 

d 

(2) the series gy{y,z) :— -—g{y,z) satisfies lim gy{y,p) = -\-oo, 

oy' y^n(p)- 

(3) the function ri{z) is continuous at p, and 

(4) the radius of convergence of q{z) — Zc{0, z"^, z^ . . .) is larger than p. 
Then the unlabelled class Q is subcritical. 

Proof. We have to show that the list of four criteria above implies that (i) p is 
non-zero, (ii) A{z) is analytic at p, and (iii) g{y,z) is analytic at {f{p),p). The first 
criterion exactly implies (i). It is actually in (0, 1) (see the remark after Definition[lO| 



about p being smaller than 1). And we have already shown in Lemma 11 that 
p e (0, 1) automatically implies that A{z) is analytic at p, which proves (ii). Next we 
show (in) holds. First we show that f{p) < r]{p), that is, y 1-^ g{y,p) is analytic at 
f{p). If ?7(p) < f{p), then g{f{z),z) is infinite a.t z = p, so f{z) = z exp{g{f{z), z)-\- 
A{z)) is also infinite at p, which is impossible (any solution of a strongly recursive 
system is finite at its radius of convergence). The case r]{p) — f{p) is excluded in 
a similar way as in Lemma [9] More precisely differentiating gives 

f'iz) = f{z)/z + (a'(z) + gy{f{z),z)f'{z) + g,(/(z), z))/(z). 

Hence, /'(z) > gy{f{z),z)f'{z)f{z) for z £ (0,p), which yields gy{f{z),z) < 
l/f{z) = 0(1) as z — p. This contradicts gy{y,p) — )■ 00 as y — )■ rj{p). Thus, 
rj{p) 7^ /(p) since the second criterion says that gy{y,z) is infinite at {ri{p),p). So 
we have f{p) < ri{p), which ensures that y i— g{y,p) is analytic at f{p). But we 
need to prove a little stronger condition, namely that g{y, z) is analytic at {f{p), p). 
Due to the continuity condition on rj{z), f{za) < ?/(z) in a small interval around 
p and therefore g{y,z) converges in a neighbourhood of {f{p),p), i.e., g{y,z) is 
analytic at (/(p), p). □ 

6. Examples of subcritical graph classes 

Recall that a block-stable graph class is completely determined by its 2-connected 
subclass B and that B must contain the link graph £ (a graph with one edge together 
with its two labelled end vertices). Let Ai — B\£. We will deal with three block- 
stable classes and show subcriticality both in the labelled and unlabelled cases: 
the class of cacti graphs, where M consists of (convex) polygons; the class of 
outerplanar graphs, where M consists of dissections of (convex) polygons; and the 
class of series-parallel graphs, where A4 consists of simple graphs obtained from 
a double edge by repeatedly choosing an edge to be doubled or to have a vertex 
inserted in its middle. Using the subcriticality criteria introduced in Sections |4] 
and [5] we will show that these three block-stable classes are subcritical and therefore 
feature a universal asymptotic behaviour with subexponential term n^^/^. 
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Theorem 15. The classes of cacti graphs, outerplanar graphs, and series-parallel 
graphs are suhcritical both in the labelled and unlabelled cases. As a consequence, 
the counting coefficient g„ of each of these classes - gn — \Q'n\/n\ in labelled case, 
Qn — \Qn\ in unlabelled case - is asymptotically of the form 

for some constants g > 0, p Cz (0, 1). The first few digits of p in labelled case and 
the approximate values of p in unlabelled case are resumed in Table\^ 

We do not claim here the originahty of the above asymptotic estimates, except for 
unlabelled series-parallel graphs: labelled outerplanar and SP graphs are treated 
in [S]; unlabelled outerplanar graphs in labelled cacti in [37] and unlabelled 
cacti in [H]). What is novel, however, is how we have derived these asymptotic 
estimates, namely through a unified (firstly for labelled and unlabelled cases and 
secondly for various graph classes) method, based on the decomposition grammar, 
the singularity analysis and the subcriticality criteria that are easy to check, which 
we will show below. Analogous asymptotic estimates hold for classes of graphs that 
are stable under taking connected, 2-connected, and 3-connected components, but 
have only finite 3-connected subclass (see Subsection 6.4). 

To prove Theorem |15[ we check whether the sufficient conditions for subcriticality 
(Lemma [9] in the labelled case. Lemma f 4 in the unlabelled case) are satisfied. 



Throughout this section we use the functions and notations in Lemmas [9] and [T4] 

6.f. Cacti graphs. The asymptotic study of cacti graphs is carried out in [37] for 

the labelled case and in [531 for the unlabelled case. Cacti graphs are such that 
M = B\£ consists of (unoriented convex) polygons. Therefore, the derived class B' 
is isomorphic to unoriented sequences of at least two vertices (because the polygon 
can be broken at the root-vertex). Thus 

B'{y) = y 



2(1 -y)' 

where the terms y counts the link-graph with a pointed (i.e. distinguished but unla- 
belled) vertex. The series g{y) = B'{y) clearly diverges at its radius of convergence 
1. Therefore the class of labelled cacti graphs is subcritical. 

In the unlabelled case, we have to take automorphisms into account. The only 
possible symmetries of unoriented sequences are the identity and the order-reversing 
of the sequence, therefore 

2(1 -Si) 2(1 -S2) 
and the series g{y, z) ~ Z^'iy, f{z^), f{z^), ■ ■ •) satisfies the expression 

From this we obtain the equation satisfied by y = fi^)- 

f{z? , l + f{z) 



y 



2(1 -/(z)) 2(1 -/(z^))- 

Since f{p) < 00 and h{z) ^ f{z) (i.e. h{z) is coefficient-wise dominated by f{z)), 
we have h{p) < 00, so that f{p) < 1. As a consequence 1/(1 — f{z^ j) is analytic at 
p. Looking at the expression of g{y,z), we see that the radius 77(2) of convergence 
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of y I— > g{y, z) satisfies 77(2) = 1 for z close to p (in particular rj{z) is continuous 
at p) and gy{y,p) goes to infinity (double pole in y) when y — >■ 1. It remains to 
check that q{z) is analytic at p. The cycle index sum of polygons is well-known 
(the automorphism group is the dihedral group for each polygon): 

ry , ^ Si 1 0(r) 1 ^ , sf + 2siS2 + S2 

ZB{si,S2,...)^-— + -}^——\0g' 



2 2 ^ r Sr 4(1 - S2) ' 

r>l ^ ^ 

where </> is the Euler totient function. Hence 

'^(-) - 4(1 -/(z^)) + 2 1^ ~ " 

This function is clearly analytic at p, which concludes the proof that the class of 
unlabelled cacti graphs is subcritical. 

6.2. Outerplanar graphs. Our second example, outerplanar graphs, has been 
studied asymptotically in ^Sj for the labelled case and in [4j for the unlabelled case. 
An outerplanar graph is a graph that can be embedded in the plane such that 
all vertices lie in the outer face. It is also defined as the class of graphs avoiding 
1^4 and K^ i as minors. A well-known characterisation of 2-connected outerplanar 
graphs with at least 3 vertices says that they are dissections of a (unoriented convex) 
polygon. So the computation shares some resemblance with the one for cacti graphs, 
except that the polygon is filled with chords in a planar way. In the labelled case, the 
classical dual construction says that dissections of oriented polygon are in bijection 
with rooted plane trees with no node of degree 2. The leaves of the tree correspond 
to the edges (minus one) of the dissection, which are themselves equinumerous with 
the (non-rooted) vertices of the dissection. Therefore 

B\y)=y+\T{y), with T(y) = , 

where the first term y in B' (y) takes account of the derived link-graph. The series 
r(y) satisfies T(y) (1 - 3y - i/y^ - 6y + l)/4, therefore it has a square-root 
singularity at its radius of convergence 3 — 2\/2. This ensures (from the 

remark after Lemma [9]) that the class of labelled outerplanar graphs is subcritical. 

Next consider the unlabelled case. We will first compute Zg/. The only possible 
symmetries are the identity and the reflection along an axis passing by the rooted 
vertex. One finds (see [3J [32] for a detailed calculation) 

^B'(S1,S2,...) = \^{si) + -^^^^(«2)- 

Hence the series 5(2/, z) — Zjs>{y, f{z'^), f{z^), ■ ■ ■) satisfies 

g[y^z)^\T{y) + y±0^T{f{z^)). 

We first observe that g{f{p,p) is finite, since it is smaller than f{p) which is 
finite, more precisely, g{f{p),p) gathers from f{p) the connected rooted graphs 
with a unique block incident to the pointed vertex. In particular T{f{p)) is finite, 
which ensures that f{p) < px- Since p < 1, /(z^) is strictly smaller than pT in 
a neighbourhood of p, so T{f{z'^)) is analytic at p. Consequently, for any z close 
to p the dominant singularity r]{z) of y H' g{y, z) is px (because of the first term 
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T{y)/2), hence 77(2) is continuous at p (it is actually constant equal to pT around p), 
so the criterion (2) of subcriticality is satisfied. Moreover, g{y, z) inherits from Tiif) 
a square- root expansion at pT, so gy(y,z) diverges as y ^ px, thus the criterion 
(3) is satisfied. 

It remains to check the criterion (4). To find an expression for one has to 
enumerate dissections of a polygon under all possible symmetries (rotation and 
reflection), where the duality with trees helps to formulate the decompositions. All 
calculations were done (see 

d>l '^2 

where P{si,S2) is a certain polynomial in si and S2. Consequently, the series 
q{z) = Zb(0,/(z2),/(2;3)^...) satisfies 

d>2 

The argument of log is non-zero around p, because for d > 2, f{z'^) < pT a. 
neighbourhood of p, and T{y) -\- y is less than 1 at its singularity pT (precisely 
it equals (5pT + l)/4 < 1). Hence the criterion (4) is satisfied. In conclusion. 
Lemma [M] implies that the class of unlabelled outerplanar graphs are subcritical 
class. 

6.3. Series-parallel graphs. We use the decomposition grammar for 2-connected 
SP-graphs developed in [8 and [T9l . More precisely, we will use the decomposition 
grammar for the so-called series-parallel networks. A network is a connected SP 
graph with two pointed (distinguished but unlabelled) vertices, which are called 
the poles and denoted by — and -|-, such that adding the edge (— , +) results in a 2- 
connected SP graph (possibly the two vertices — and -I- are already adjacent in the 
network). We call series-networks those that are not 2-connected (when deleting 
the edge between the poles if any) and parallel networks those that are 2-connected 
with at least 2 edges. The classes of networks, series-networks and parallel-networks 
are denoted respectively V, S, and V. The link-graph consisting of one edge from 
— to -I- is denoted by e to distinguish it from the link-graph i consisting of an edge 
and two labelled end vertices. 

In the labelled case we obtain the decomposition grammar and its corresponding 
system of EGFs (on the right is the associated system for the EGFs, Set>fc means 
Set constrained to have at least k components): 

V = e + S + V, r V{y) = l+Siy)+V{y), 
S = V-Z-{e + V), Siy) = yV{y){l + V{y)), 

V = e-Set>i(5) + Set>2(5) [ V{y) = 2exp(5(y)) - 5(y) - 2 

The system for the EGFs is clearly strongly recursive and the series I?(y) is 
easily aperiodic. In addition, the function-system r(yi, ?/2, J/s, y) with j/i = 'D{y), 
2/2 = S{y), 2/3 = 'P(y) defined by the right-hand side of the system is clearly 
analytic everywhere (because exp is analytic everywhere). In addition, easy lower 
and upper bounds imply that the radius of convergence ryu of 2?(y) is in (0,1). 
Hence Theorem [2] applies, ensuring that all the series for networks have a square- 
root expansion at their common radius of convergence rjx>- Now we want to show 
that g{y) := B'{y) has radius of convergence 7725 and that g'{y) goes to 00 as ?/ — > r;^. 
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First, note that in the labelled framework the number of networks is at least as large 
as the number of vertex-rooted 2-connected SP graphs (because each graph in B' 
gives rise to c? > 1 networks, with d the degree of the root- vertex) . So the series 
y^'D{y) dominates coefficient-wise the series yB'{y), shortly written B'{y) ^ u'D{y). 
Surprisingly the domination also goes the other way. Precisely speaking, we have 
the inclusion -V (-B' , as illustrated in Figure [l] so that y^'D{y) ^ B'{y). 




Figure 1. Injection from networks to vertex-rooted 2-connected 
SP graphs by adding one vertex and two edges. 



We thus have y^'D{y) ^ B'{y) ^ u'D{y), which ensures that B'{y) has the radius 
of convergence rjj} and that g{y) = B'{y) satisfies g'{y) — >■ oo as y — 77^, since ^{y) 
has a square-root expansion and hence f (y) — >■ cxd as ?/ — > 77^). This concludes 
that the class of labelled SP graphs is subcritical. 

Now we turn to the unlabelled case, which is technically more involved. For 
each class J\f £ {VjS./P} of networks, the automorphisms have to fix each of 
the two poles + and — , and the cycle index sum Z/v(si, S2, ss, • ■ •) is defined as 
a sum of monomials over all automorphisms, as in Section 2.1 We define N = 
N{y,z) :— Z^{y, f{z'^), f{z^), . . .). For networks, the automorphisms have to fix 
each of the two poles -I- and — . However, it is also useful to consider networks up 
to exchanging the poles and their corresponding series Za/'(si, S2, S3, ■ • ■) defined 
similarly as the cycle index sum, but summing over the automorphisms exchanging 
-|- and — (instead of those fixing -I- and — as in the series Z^). Define the bivariate 
series N = N{y,z) := 'Zj^{y,f{z'^),f{z^),...) for each class N € {'D,S,r} of 
networks. For a bivariate series f{y,z) and fc > 1 define f''^\y,z) :— f{y'',z''). 
The decomposition for pole-fixing networks and pole-exchanging networks give rise 
to explicit systems for the cycle index sums. Under the specialization (si = y,S2 = 
f{z'^),S3 = /(z"^),...), these yield the two systems (where the argument {y,z) of 
the series is omitted): 



(PF)<^ 



( D = l + S + P, 
S = yD-{l + P), 
P = 2exp(5:^^^V))-5-2, 



and for the pole-exchanging networks 

D = l + S + P, 
S = Di^) ■{y + f{z^){l + P)), 



-S 



(2fc-l) 



2k 



2k -I 



S-2. 
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The series g{y, z) Zjs' {y, f{z'^), f{z^), ■ ■ •) can be expressed in terms of the series 
in the above two systems (this is done in |19j using the dissymmetry theorem): 

giy,z) = ^y\l + P)D+^f{z^)D^^^-{l + y + yP^^^+p) + 

2y exp (y2^^, - '^v - ^vS - ^y^' - ^y^^'^ - ySP. 



As in Lemma 



14 



we let f{z) ■.= C'{z), g{y,z) := ZB,{yJ{z^)J{z^),...), and let p 
be the radius of convergence of f{z) and r = f{p). In addition, we let r]o{z) and 
r]jy{z) be the singularity functions of y i— >■ D{y, z) and y t-^ D{y, z). 

Lemma 16. Each of the series D{s, z), S{s, z),P{s, z) admits a square-root expan- 
sion at {riDip),p)- 

Proof. As Zigi{f{z),f{z'^), . . .) counts derived connected SP graphs with a unique 
block incident to the pointed vertex and f{z) counts all rooted SP graphs, we have 
9if{p),p) = ZB'{f{p)J{p^),fip^),...) < f{p). In addition f{p) < oo since f{z) is 
a solution of nonlinear functional equations. Hence g{f{p),p) < oo. As ■ V C B' 
(see Figure [T]) we have 

u^D{y,z) di g{y,z). 

Hence, by domination, D(f{p), p) < oo and therefore the singularity of y M- D(y, p) 
is larger than /(p), which implies in particular that z = p is an admissible valua- 
tion. The equation-system (PF) F(y; z) for y = {D,S,P,y) satisfies all conditions 
required in Theorem |2] (aperiodicity of D, strongly recursive system, analyticity of 
the functional system F(y; z) at the required point). Therefore the series D, S, and 
P admit singular expansions at {r]r){p), p). O 

Lemma 17. If rjulp) < rfjj(p) then the series D{s, z), S{s, z), P(s, z) are analytic 
at {rjoip), p) , if not^ then these series admit a square-root expansion at {rj-jj^p) , p) . 

Proof. Note that Figure [T] yields an injection from the pole-exchanging automor- 
phisms of T) to the automorphisms of B' . As a consequence, S2^d ^ ^B' so that 

z'^D{y,y) d g{y,z). 

Hence, Zs'{f{p),f{p'^),f{p^),...) — t < oo implies that D{f{p),p) < oo. There- 
fore r]jj{p) > f{p) so that z = p is an admissible valuation for D{y,z). Note that 
the functional system (PE) G(y; z) with y = (Z?, S*, P, y) is clearly analytic at a 
given point (Z?, S, P, y; z) (where Z?, S, P are seen as independent variables) if and 
only if D is analytic at (?/^,z^) and / is analytic at z^. Since p < 1 and rixi{z) 
is decreasing in z, G is analytic at any point such that \y\ < rii:i{p) and |z| < p. 
Therefore, the only cause of singularity for D in this domain is a branch point (i.e., 
a solution of the singularity system). From Theorem [2| we conclude that in such a 
situation, D, S, and P have a square-root expansion at {r}j=j{p) , p) . □ 

Lemma 18. The class of unlabelled series-parallel graphs is subcritical. 

Proof. We need to check that all four criteria in Lemma [M] are satisfied. Easy 
upper and lower bounds imply that p e (0,1), so the criterion (1) is satisfied. 
Let ri{z) the singularity function of y i— )■ g{y,z). Expression for g(jj,z) in terms 



■^This case is very unlikely to happen, but discarding it would require some numeric computa- 
tion of the functions t]£) and rjjy. 
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of the series of networks (pole-fixing and pole-exchanging), ensures that i](z) — 
min{rijj[z),rfj=)(z)} for z > 0. Hence ?7(z), as the minimum of two continuous 
functions, is also continuous at p, so the criterion (3) is verified. If ri{p) = rjoip) (the 
most likely case), then limj^^^^p)- gy{y,p) = oo because y'^D{y,z) ^ g{y,z) and 
because the square-root expansion of D(y,z) at {rj{p),p) yields \miDy{y, p) — oo 



as y — ^ ri{p). If r]{p) = rfjj{p), then applying Lemma 17 D{y,z) has a square- 
root expansion at {ri{p),p). Again we have Irm gy{y,p) = oo as y — ^ ri{p) because 
z'^D{y,z) :< g{y,z) and because the square-root expansion of D{y,z) at {ri{p),p) 
yields \vaiDy{y, p) = oo as y — > rj{p). Thus the criterion (2) is verified. 

It remains to check the criterion (4) . From the work in [H [19] one can extract an 
expression for K{y, z) := Z^iy, f{z^),f{z^), . . .) (longer than the one of g{y, z), but 
of the same aspect) in terms of the series {D, S, P, D, S, P}. This ensures that the 
singularity function r]H{z) of y i— >■ H{y,z) satisfies rjH{z) = min{77£)(z), 77jj(z)} = 
ri{z). Since rfnip) > f{p) and since rjn is continuous at p, we have r]H{z) > 
for z around p. Hence -ff (0, z) is analytic at p. In addition one can show that 
Zc{0, z^, z^, ■ ■ ■) — Zb{0, f{z^),f{z^), . . .), see (a combinatorial interpretation 
is that if there is no fixed vertex in an automorphism cr of a connected graph, then 
there is a unique 2-connected block B fixed by a, and the automorphism induced 
on B has no fixed vertex as well). Hence q{z) — Q(0,z) = H(0,z), which ensures 
that q{z) is analytic at p. □ 



6.4. Graphs defined by their 3-connected components. The same arguments 
used in the enumeration of the unlabelled SP-graphs can be adapted to get similar 
results for families of graphs which are defined by a finite set of 3-connected com- 
ponents (see [21] for a proper definition of these families). In these families, the 
strategy used is the same as in SP-graphs: we study the enumeration of networks in 
order to get the asymptotic counting of the graphs with a lower level of connectivity. 
The main difference is that there is an additional line in the equation-system for 
networks, which involves the Walsh polynomial (roughly speaking, the cycle index 
sum associated to the finite family of 3-connected components; see [TH] [H] for a 
detailed definition). The Walsh polynomial has two kinds of variables: ei, . . . , 
and cTj", . . . , e^r. In the equation-system for networks, these variables are to be sub- 
stituted by D^^\ . . . , Z?^'') and by -D^^"*, . . . , D^''\ respectively (in SP-graphs this 
polynomial vanishes, because SP-graphs do not have 3-connected components). 
Since the Walsh polynomial is an entire function in all its variables and has posi- 
tive coefficients, the system remains positive and the argument of Lemmas |17[ |18| 
can be adapted to this general setting in order to assure a square-root singularity 
for the counting series associated to networks, and the rest of the results can be 
easily adapted to this general framework. 



7. Limit laws 

In this section we study the limit distribution of a graph parameter defined on 
a uniformly distributed random graph on n vertices. More concretely, let Xn be 
the associated random variable. We show that all the parameters under study 
invariably converges in distribution to a normal limit distribution in the form 
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where EX„ = + 0{l) and VarX„ = a^n + 0{1) for computable constants 
fi, (J^ > 0. The main results used in this section are Theorem 2.23 of [TT] (adapted 
in this paper as Theorem |3| . In order to decide whether > we can either 
calculate it directly or use Lemma [4j For example, we can check the conditions 
of Lemma |4] with the help of small values of n, m, k. This applies for all kinds 
of graphs and parameters we consider here. However, in any case (if cr^ = or 

> 0.) these parameters are concentrated around its expected value. 

The parameters discussed in this section are the number of edges, the number 
of blocks and the number of cut-vertices. All these parameters can be encoded in 
the generating function framework of the form 

g{z,v)^exp{C{z,v)), C'(z,v)^F{C'{z,v);z,v) 
for the labelled case, and 

g{z, v) = exp {Y.k>i C{z\v'')/k) , C'{z, v) = F[C'{z, v)-z, v) 

for the unlabelled case. In both cases, z marks vertices and v marks the parameter 
under study. The singularity curve z = p{v) for the GF associated to general graphs 
is the same as the one for connected graphs, because both functions exp and exp 
are analytic functions. As a consequence, the same limit law holds for both general 
and connected graphs (hence we can restrict our study to connected graphs). In 
order to study these limit laws we analyse the characteristic system 

dF 

(15) y^F{y;z,v), l^—-{y;z,v). 

oy 

Observe that we recover the univariate system y — F{y;z), 1 — dF/dy{y; z) by 



evaluating (15) at v = 1. For the classes under study we have shown that this 
univariate system of equations is analytically well-founded (in the sense of Defini- 
tion [ij, hence by Theorem [2] y{z, 1) is the unique solution and has a square-root 
expansion around its smallest real singularity. This expansion can be extended to 
a real neighbourhood of u = 1 with expression y{z,v). 

In this section we obtain the exact parameters for the labelled case. Unfortu- 
nately in the unlabelled framework there are no explicit formulas for the expectation 
and the variance, hence we content ourselves with the normal condition derived from 
Theorem [3] combined with Lemma |4] Finally, local limit laws are also discussed for 
these parameters. 

7.1. The labelled case. The labelled case have been partially treated in f2T|. In 
this section we rediscover these results using a different and more general argument. 
Recall that t = C*{p) is the unique solution of the equation yB"{y) = 1 on the 
region where B(jj) is analytic, and p is the radius of convergence of C{z). 

We use a secondary variable v to mark edges. For a graph class G we write 
G = ^n,mGn.rrn whcrc 71 deuotcs the number of vertices and m denotes the number 
of edges. We define the bivariate generating function as 

G{z,v)=y^-\Gn,m\z"v"'. 

The series is exponential according to the (labelled) vertices and ordinary according 
to the (unlabelled) edges. A similar definition can be stated for both connected and 
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2-connected graphs in Q (C(z, v) and B{z, v) respectively). For a block-stable graph 
class Q, the block-decomposition translates to 

C'{z,v) = z exp{B'{C'{z,v),v)), 

where B'{z,v) denotes the derivative of B{z,v) with respect to the first variable. 
Observe that we need to deal with the bivariate GF associated to 2-connected graphs 
of the class. All subcritical classes under study are subclasses of planar graphs, 
hence the number of edges is linear. Consequently the subcritical condition is stable 
under a slight variation of v around v = 1. We denote by B'^{z, v) = ^B'{z, v), the 
derivative of B'{z,v) with respect to the second variable. 

The next parameter considered is the number of blocks, which is coded by the 
following equation 

C'{z,u) = z exp {uB'{C'{z,u))) , 

where the secondary variable u marks the number of blocks. Finally, for the number 
of separating vertices, also called cut-vertices, we consider the bivariate GF 

C'{z, w) = zw (exp {B' (C*(z, w))) - 1) -t- z. 

Observe however that w does not mark exactly the number of separating vertices, 
as the root-vertex is always considered as a cut-vertex. However, the limit law for 
the number of cut-vertices does not depend on the behaviour of a single vertex. 

All the results are contained in Table [2] Despite the methodology used here 
(direct application of Theorem [S]) is different from the one of [51] , the parameters 
we obtain are the same. 



Parameter 


F{y;z,v) 






Edges 
Blocks 
Cut-vertices 


zexp{B'{y,v)) 
zexpivB'iy)) 
zv (exp(S'(y)) -l) + z 


log (^) 


log 


-, (re;'(r,i))^ 

^ 1+t^B"'(t) 
(p) l+r2B"'(r) 
p ^ 1+t^B"'{t)) 



Table 2. Parameters with the corresponding value of ^ and . 



7.2. The unlabelled case. The same results used in the labelled framework work 
for GFs in the unlabelled case, and Theorem [3] assures the asymptotic normal limit 
distributions (with a non-vanishing variance) for these parameters. However, we 
cannot obtain closed expressions (as the ones in Table [2]) for both the expectation 
and the variance of the resulting random variables, due to the Polya operator. 

Using the techniques introduced in Section [9] one can find numerical approxima- 
tions for the derivatives of the function F{y] z, v). Such computations are conducted 
in [4] for the number of edges in random unlabelled outerplanar graphs. 

7.3. Local limit laws. In addition to the central limit laws described above, (much 
stronger) local limit laws hold for the parameters under study, either in the labelled 
or the unlabelled framework. The main results used here are ^ and the asymptotic 
estimates for both [z"]C*(z) and [z"]C*(z). For conciseness, we state the result for 
connected graphs, but a similar result is also valid for general graphs. 
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Theorem 19. Let Q be a subcritical graph class with the property that [z"]C{z) > 
for n > Hq. Let X.^ be a random variable for a graph parameter studied above ( the 
number of edges, blocks, or cut-vertices) defined on a random uniformly distributed 
connected graph of Q on n vertices and assume that > 0. Let m = ^n-\-0{n^/^). 
Then 

8. Degree distribution 

In this section we discuss the degree distribution of graphs of a subcritical graph 
family. We denote by and X^ the random variables which count the number of 
vertices of degree fc in a randomly chosen member of size n of the family Q and 
respectively. We further denote by dfe and the probability that the root vertex 
of a randomly chosen member of Q' or Q' is k. The main tool to obtain asymptotic 
results is the analogue on systems of equations of Theorem [3] 



.1. The labelled case. Consider the series 



n! 



where we use the notation v = (wi, . . . ,wa;,Woo) and &j.„.„j „^ „^ is the number 
of derived 2-connected graphs with 1 + n — 1 + ni + ■ ■ ■ + Uk + n^o vertices such 
that one vertex of degree j is marked and the remaining n vertices are labelled by 
1, 2, . . . , n and ng vertices have degree ^, 1 < < fc, and vertices have degree 
greater than fc. By definition we have 

B'{z)^ B'^iz.l). 

2<j<oo 

Hence the radius of convergence of the functions Bj{z, 1) is greater or equal than 
the radius of convergence of B'{z). 
We introduce an analogous series 

C'j(2,v)= S;";"l,---."(:,"oo''^l ^ ' ' ' ^fc'°'^oo°° ^) 

n,ni , ... ,nfc 

for derived connected graphs. For convenience, we set Cq(z,v) = 1, which cor- 
responds to the case of a graph consisting of a single derived vertex. We further 
set 

B'{z, w)^Y. ^'A^' 1)^' = E 

and 6^ ^ is the number of derived 2-connected graphs with n + 1 vertices such 
that one vertex of degree j is marked and the remaining n vertices are labelled by 
1, 2, . . . , n. Analogously, we define the function C'{z, w) = J2j ^'ji^^ Accord- 
ing to the block decomposition of connected graphs 

C'{z,w) = exp(B'(zC"(z),u))). 

In the following, we set Bj{z) := Bj{z, 1). We will also use the series B(z, v), which 
is the refined version of the generating function B(z) of unrooted blocks taking into 
account all vertex degrees. 
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Theorem 20. LetQ he a family of random subcritical graphs (such that [z^^]C{z) > 
for n > Uq) and Q' be its derived family. Then for k fixed, 

(1) the limiting probability dk that the pointed vertex of a member of Q' has 
degree k exists and is given by 



(2) the random variable that counts the number of vertices of degree k in a 
randomly chosen member of Q satisfies a central limit theorem with mean 
EX,^' = dfcft + 0(1) and variance VarX^ = o-|.n + 0(1). 

Proof of Theorem^^. We use the following lemma, which makes it easy to com- 
pute dk. The proof is the same as the one for [TI] Theorem 9.17]) 

Lemma 21. The generating function p{w) ~ "Y^dkUi^ satisfies 
p{w)^pe^'^^-^'>^B'{z,w) 

z=pC'(p) 

and p{l) = 1, thus the dk 's define indeed a probability distribution. 
Therefore, we have 



4 = [w'^Mw) - [i«'=]pe^'(--)|-S'(z,i«) 

oz 



z=pC'{p) 



z=pC'(p) 



\i=l V k 



,k\2 



z=pC'(p) 



Ea7^^'W|E E B,Az)---Bi„Az) 



z=pC'[p) 



From there the result follows, as the second term is exactly the representation of a 
connected graph of root degree k—i according to the block decomposition, evaluated 
at z = p. □ 

Proof of Theorem^^. First we derive a system of functional equations which is 
satisfied by = Cj (z,v), and is refinement of 

C*(z) = zexp(B'(C*(z))). 

Lemma 22. Let Wj = Wj{z,v;C[, . . . ,C[,C^), j e {1,2, ... oo} defined by 

k 



k-j 
i=0 



CKz,v) + C^(z,v) 

Ki=k-j + l 



l<j<k, 



W^=v^ K]c:(z,v) + C^,(z,v) 



\i=0 
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Set W = {Wi, . . . , Wk, Woo)- Then the series C^, . . . , C^, satisfy the system of 
equations 

c'^ (z, v) = exp ("^ i?; (z, w) + s;, (z, w) V 1 - n ^^^^^ . 

\j = l / l<£i+---/£ffc<fc r=l 

Proof. As already indicated, the proof is a refined version of the functional equation 
fulfilled by C*(z), which reflects the decomposition of a derived connected graph 
into a finite set of derived 2-connected graphs, where every vertex (different from 
the root) is substituted by a derived connected graph. The functions Wj serve 
the purpose of marking (recursively) the degree of the vertices in the 2-connected 
blocks which are substituted by other graphs. In the definition of Wj, the summa- 
tion means that we are substituting a vertex of degree i, but since originally the 
vertex had degree j, we are creating a new vertex of degree i + j, which is marked 
accordingly by Wi+j. The analogous holds for Woo- O 

To prove Theorem [20] ([2| we observe 

0<j<oo 

Furthermore, since the above system of equations is strongly connected, all func- 
tions Cj{z, 1) have the same radius of convergence as C'{z), as mentioned in Section 



3.2 By assumption this radius of convergence is smaller than the radius of con- 
vergence of B'{z). Hence, by stability if v is sufficiently close to 1 then the the 
singularities of Bj and C'j do not interfere, in particular we can apply Theorem [2] 
and obtain that all functions have a square-root singularity. Finally, let 

n! 

n.m 

be the generating function for the numbers Ck-^n,m of unrooted connected outerpla- 
nar graphs of size n with m nodes of degree fc. Then C'-'^-'(z, v) satisfies 

g!^'^^ = E 1' ■ • ■ ' 1' ■f'' 1) + v(^kiz, 1, • • • , 1, f , 1) + C'ooiz, 1, . . . , 1, 1), 

and thus C'^''^(z, v) has a singular expansion of order | around v — 1- Furthermore, 
the central limit theorem for the number of vertices of given degrees with mean /i^n 
and variance a^n follows by the analogue of Theorem [s] for systems of equations. 
It immediately follows that dk = fJ-k as there are exactly n possible ways to root 
an unrooted object of size n at one of the vertices and thus the probability that 
a random vertex has degree k is exactly the same as the probability that the root 
vertex has degree k. Despite that, we could also use formula ([t]) to compute /x^, we 
would obtain the same value as for dk here. 

□ 
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8.2. The unlabelled case. We introduce cycle index sums 

Zg,Xsi,S2, Wl,2, ;VkA,Vk,2, ■ ■ ■] Woo,l J "^00,2 , • ■ •) 

for the class of pointed blocks, where the pointed vertex has degree j and is not 
counted, and where the variables Vij count the cycles of length j of vertices of degree 



i, and Woo.j counts those vertices of degree greater than k. As in Section 8.1 let v = 
{vi, . . . ,Vk,Voo)- Denote the corresponding OGFs by Bj{z,-v),j G {!,..., fc, 00} 
and let 

k 

B\z, v) J2 ^) + ^ooB'^i^, v). 

Note that 

Zg,{si, S2, . . . ; 1, 1, . . . ; . . . ; 1, 1, . . . ; 1, 1, . . .) = Zg,{si,S2, ■ ■ .), 

and thus the singularity pi(v) of B'{z, v) is the same as that of B'{z) at v = 1, and 

/9i(v) is the dominant singularity of the system B'(z,v) = {Bj{z,-v))jfz[i k,oo}- 

We now introduce the multivariate generating functions 

n;ni,...,nfc,noo 

where the coefficient Ci.n;nx,...,nk,n^ denotes the number of elements of size n of C, 
where the pointed vertex has degree j and with rij, i e {1, . . . , fc}, vertices of degree 
i and rioo vertices of degree greater than k. We further set 

B'{z,w) =Y.B',{z, l)w^ = E^",.^"^' 

and 

j n,j 

As we need cycle indices for the block decomposition, we set 

Z^,{S1,S2, . . . ; W) = E (■51: S2, ■ • O^^'^- 

i 

Note that the variable w, which counts the degree of the root, is not involved in 
any permutation cycle. Then, 



C'(z,zi;)=exp lj2]zg,{z'C'{z'), 
\e>i 



z^'C'{z^'),...-w') 



Theorem 23. Let Q he a family of random subcritical graphs (such that [z"] C{z) > 
for n >nQ) and G' be it's derived family. Further, let dk be the limiting probability 
that the root vertex of a member of Q' has degree k and let be the random variable 
that counts the number of vertices of degree k in a randomly chosen member of Q . 
Then 

(1) dk = p (E-=i ^^Z^,{z.p'C'{p%p^C'{p% . . ■)\^^^cyf'k-^ip)) ■ 

(2) satisfies a central limit theorem with mean = pt^n + 0(1) and 
variance Var X'^ — tr^n + 0(1). 
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Remark 24. In the unlahelled case, we cannot expect that jik equals dk, as vertex- 
rooting is only possible at fixed points of permutations, and thus in unlahelled graphs 

nth. + g'n- 

Proof of Theorem ([T]) . The proof is based on the following lemma: 
Lemma 25. The generating function p{w) = ^ d^w^ is equal to 



p— exp I Z^,{u,z^C'{z^),...;w) + 
on 



£>2 ) 



and p{l) = 1. That is, it defines indeed a probability distribution. 
Proof. We use til] Lemma 2.26] with 

f{z,w) = zC'{z), 



Y,\z^,{/C'{z%...;w')\ 



H{z, w, u) = exp ^^g, (u, z C'{z ),...; w) + 

and the same line of reasoning as in [TTJ Theorem 9.17] proves the first part of the 
lemma. For p(l) = Hu{p, l,pC'{p)) we obtain 



pC\p)^Z^,{u,p^C\p%...) 



= 1 

u=pC'(p) 



by the implicit function representation of C'{z). □ 
Now we can determine d^ = [w'']p{w), which is equal to 



e^P ( E UB'ip'C\/), ■.■■,w')] (^Z^,{u, p'C'ip'), ...;w) 



Recalling Zg, = X]/>i ^b' ' obtain that the previous expression can be written 



^B' ^ ^1>1 

P i2^ZB'S^'P"^'(p")'---)\u^pC'ip) ■ [-'"^]^-P (^^]zg,i/&{/),...;w')^ . 



Observe that the second term translates into C'^_^{p), as the following equality 
holds: 



e^p\J2^^Zg,{p'C{/ 
\e>i * 



□ 



Proof of Theorem As in the labelled case, we observe that the functions 

C'j{z,v) satisfy a system of equations, using a refinement of the block decomposi- 
tion. In the following we will denote by J2i^ — J2i=r + 
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Lemma 26. For each j — 1,2, . . . ,k, oo, let Wj be defined by 

k-e / k,oo 

z^O \i^k~j+l 

Yl^C[{z,x)\ . 
Set Wj,i = W-j{z\ vl,vl^) and 

W« = iWi,i, Wi,2i, Wkj,Wk,2i, . . . ; W^,i, W^,2i, ■ ■ •)• 

We denote by 5„(si, S2, . . .) the cycle index of the symmetric group on n elements. 
Furthermore S'„[Zb] denotes the substitution s; ^ Ze(s;, S2;, . . . ; W*-'-*), Z > 1 in 

S„{si,S2, ■ ■ ■)■ _ _ 

Then the series C[, . . . ,C'f., C'^ satisfy the system of equations 



ll+2l2 + ---+jlj=J r=l 



I , 

J>1 \r=l 

j 



l<li-\----\-klt,<k r^l 



, ^ , . . . , 



Remark 27. T/ie functions Wj^i and thus the whole system can also be considered 
in terms of cycle index sums, using the cycle index sum 



l'fc,l,Wfe,2,---;Woo,l'^oo,2,---) 



for rooted connected graphs. The root vertices in Wjj are fixed and thus have cycle 
length 1 . We will need this terminology in the proof of Lemma \29[ 

Proof. As in the labelled case, we refine the recursive decomposition of graphs into 
it's 2-connected components. The functions Wj^i plays the analogous role as in the 
labelled case, except that we need a second index for representing the cycles of 
different length appearing in the cycle indices. This directly leads to Cj(z,v) for 
j = 1, . . . k. For C'^{z, v) we obtain 



E n^'- 

ih,h, ■■■,lk,l^) r=l 



E Y[S,\z,,(z,z^,...;W('^ 



h+...+klk<k r=l 



where the first sum rewrites into the exponential term appearing in C^(z,v). □ 
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It is easily checked that the system is strongly connected as every Cj{z,v) 

depends on C^(2,v) for j e {!,..., fc} and C'^{z,v) depends on all the values 
j £ {1, . . . , fc}. Obviously, 



(16) C'(z) = 

Define 



k 



v=l 



C'(z,v) :=^«,C;;.(z,v)+«ooC^„(z,v). 

3 = 1 

Then the generating function of derived connected graphs of Q, where the variable 

V counts the nodes of degree k, is given by C"(*'^(z,v) — C'{z,Wk), where — 

Liemma 28. C'^^^{z,v) has a square-root singular expansion around its singularity 
P2{v) in a neighbourhood of v = 1. 

Proof. By Equation ([l6| the singularity of the system C'(2;,v) at v = 1 is p2, 
the same as that of C'{z). As the system is strongly connected, every Cj(2:, 1) 
has radius of convergence p2- Since C'''^\z,l) is a linear combination of these 
functions, it has the same singularity, which fulfills p2C'{p2) < Pi due to the sub- 
criticality assumption. By the stability property of subcriticality it follows that 
((7'W(z,u),i?'('=)(z,w)) is subcritical near 1, and hence we obtain a square-root 
singular expansion. □ 

Consider the cycles index sums 

%(S1, S2, . . . ; Wl,l,'Ul,2, ; Wfe.l, Wfe,2, • • • ; Wooa, Woo,2, • • •) 

for (unpointed) connected graphs of Q. By taking si — s, Si ^ z^ for i > 2, 
Vj^i = 1 for I < j < k,i > 1 and Vk.i — for i > 1 we obtain its corresponding 
OGF (7('^'(z, v) = Z^{z, z^, . . . ; 1, 1, . . . ; . . . ; w^, . . . ; 1, 1, . . .), where the variable 

V counts the nodes of degree k. 

Liemma 29. C'-'^^z^v) has a singular expansion of order | around its singularity 
P2{v) in a neighbourhood of v ~\. 

Proof. We have to express the system of equations in Lemma |26] in terms of cycle 
index sums and analyze the trivariate generating functions 

C-e^) (s, z, i;) = %(s, z2, z^^ . . . ; 1, 1, . . . ; . . . ; v, . . . ; 1, 1, . . .). 

Obviously, C^^\z,z,v) = C^'^^z.v). Analogously we define C'^^\s,z,v). We ob- 
tain 

C('=)(s,z,«) = C('=)(0,z,w) / C'('=)(s,z,i;)ds, 



by the same arguments as in Section [sj Due to stability, we obtain a square- 
root singular expansion for (7'^'^-'(s, z, w) as before, with a singular term of the 
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form (1 — s/p{z,v)y^^. Integration and subcriticality conditions lead to a singular 
expansion of the form 

C''''\s,z,v) = a{s,z,v) + 6(s, z,v) (l - — — - ) 
At s = z we can represent the singular part as 

3/2 / ^ \ 3/2 

k{z,v) ( 1 



with an analytic factor k,[z,v) and some analytic function p{v). Since the singular 
manifold of C'^*-') (s, z, v) is the same as that of C^^'^ (z, v), that is z — p{z, v) if and 
only if z = /02(i'), it follows that p{v) — P2{v)- D 

Now, with the singular expansion in p{v) given, we can immediately deduce a 
central limit law by the analogue to Theorem [3] for systems of equations. □ 

To calculate flk we can use Expression ([?]). The derivatives give 



/c,oo k 
3=1 i=l 

bTF,(C(p,l);p,l)= \E-Y,{B'^),{z,W) 



(C(p4);p4) 



(C(p4);p4) 



where yu-i denotes the (fc — i)-th coordinate of y satisfying y = F(y;z,v), {B'j)y. 
is the derivative of B^ with respect to Vi and E is the exponential term appearing 
in C'^. As the formula includes all partial derivatives of (Bj), we see that the 
calculation of pk will be very involved and it is (in general) not equal to dk- 

9. Computing the growth rate 

In this section we show how to compute the exponential growth of unlabelled 
2-connected, connected and general SP-graphs, respectively. The main tool is the 
decomposition grammar and the systems of functional equations developed in Sec- 
tion [6|3l 

We will first review the growth constants for subcritical classes of graphs in 
Section |9.1[ Next we will approximate the exponential growth for the number 
of unlabelled SP-graphs in two steps: the first step for 2-connected unlabelled 



SP-graphs (Section 9.2 1 and the second step for unlabelled connected SP-graphs 
(Section 9.3 1. To this end, we start with a functional system of the form y = 
FAr(y;z), which is a truncated version (combined with iteration) of the original 
functional system F(y; z) determining y and provides an approximation of y. The 
solution of this system, together with the additional restriction given by a singular 
equation (given by the determinant of the Jacobian of the previous system) gives 
the desired growth constant. 

The error introduced in the computations becomes smaller as soon as we take 
a better approximation for y (in other words, more terms on the Taylor series of 
its components). We refer the reader to [1] and [5S], in which the authors study a 
similar iterative scheme. 
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9.1. Explicit growth constants. The computation of the growth constant for 
the classes under study is based on the subcritical condition. In labelled case, the 
singularity p is given by the equation p = r exp (— y8'(r)), where r is the small- 
est solution of the equation yB"{y) — 1. In case of unlabelled case, the growth 
constants obtained earlier in the literature are based on explicit expressions for 
Ze'isii S2, ■ ■ ■ ), e.g. the classes of cacti graphs and outerplanar graphs. Here we 
will show how to derive the growth constants of SP graphs, even when there is not 
an explicit expression for Ze'(si, S2, ■ ■ ■)■ 

We have shown in Theorem 15 that the asymptotic number of cacti, outerplanar 
and SP graphs is of the form c n^^/^ n\ (1 + o(l)) in the labelled case and 
is of the form cn"^^^ p^" (1 + o(l)) in the unlabelled case. In both cases (either 
labelled or unlabelled) the subexponential term n~^/^ suggests that all these classes 
have an arborescent structure. In Table [s] the exponential growth constant p~^ for 
connected (and general) acyclic, cacti, outerplanar and SP-graphs are shown. 



Family 


Labelled 


Unlabelled 


Acyclic 


2.71828 


2.95577 


Cacti 


4.18865 


4.50144 


Outerplanar 


7.32708 


7.50360 


Series-Parallel 


9.07359 


9.38527 



Table 3. Exponential growth for distinct subcritical classes. All 
constants are referred to connected and general classes. 



9.2. Parameters for 2-connected SP-graphs. First we compute the growth 
constant using the equations for networks (and related classes). Later, the re- 
sulting singularity is transferred to the counting series of 2-connected SP-graphs 

ZBiz,z'^,Z^,...). 

Denote by Zs{si, S2, ■ ■ ■ ), Z-p{si, S2, ■ ■ ■) and Zx>{si, S2, . . . ) the cycle index sums 
associated to series, parallel and general networks. Let S{z) = Zs{z, z'^ , . . .), 
P{z) = Z-p{z, z'^, . . .), D{z) = Zx)(z, z^, . . . ). As it is shown in Section [6j the 
system of equations which defines series, parallel and general networks is 

r D{z) = 1 + S{z) + P{z), 
(17) \ S{z) = zD{z){l + P{z)), 

[ P{z) - 2exp(E.>a^«W) -5(z)-2. 

Let U{z) = exp (X)i>2 ^S'^^''{z)]. This function is analytic at the singularity of 



S{z). We can isolate S{z) from System (17 1 and obtain the implicit relation 



S{z) = z (2exp(S'(z))C/(z) - 1) (2 exp(5(z))C/(z) - 1 - 5(z)) . 

Observe that U{z) depends on S{z). In order to get an approximation of the 
smallest singularity of S{z), we start looking for an approximation of S{z). This 
can be done in the following way: let So{z) — and Uo{z) — 1. Let be a 
positive integer. For a function /(z), analytic at the origin, define pol^ {/(z)} as 
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the Taylor polynomial of degree N of f{z). Define Sk+i{z) and Uk+i{z) recursively 
in the following way: 

Sk+i{z) = polj^{z{2eMSk{z))Ukiz) - 1) (2 exp(5fc(z))C/fc(z) - 1 - Skiz))} , 
U,+,{z) = pol^{exp(E.^a4+iW)}- 
We stop when Sk~i{z) — Sk{z) (equivalently, when Sk{z) = pol^{S'(z)}). Denote 
then u{z) — exp (^J2iL2 l^k\-^)^^ consider the solution s{z) of the equation 
s{z) — z (2 exp(s(z))w(z) — 1) (2 exp(s(2))it(z) — 1 — s(z)) . 

Taking u(z) instead of V(^z") introduces an error, which can be controlled in that the 
larger iV is, the smaller the error. Under these assumptions, s(z) is defined by an 
equation of the form s = H(s;z), with H{y; z) = z(2 exp(y)u(z) — 1)(2 exp(y)u(z) — 
1 — y). Consequently, we find its smallest singularity by solving the characteristic 
system y — H{y; z), 1 — Hy{y; z). In Table|4]the values obtained for several choices 
are shown. In particular, one can observe that the accuracy in the computation are 
improved by increasing the order of truncation. 



N 


Singular point {p) 


5 


0.12421863192426192376 


10 


O.12419919715484630978 


20 


0.12419909378526277564 


50 


0.12419990937841528588 



Table 4. Values for N and the corresponding value for p. 



These computations give the truncated value of z, pi w 0.12420. This constant 
is slightly smaller than the one which is obtained in the labelled case (whose value 
is approximately 0.12800. See [5]). This singularity is the same for the parallel 
family and for general networks. Using the same ideas one finds that the smallest 
singularity for exchanging poles networks is strictly bigger than pi. Consequently, 
applying the transfer theorems of singularity analysis we conclude that the number 
hn of 2-connected unlabelled series-parallel graphs on n vertices is 

where 71 — p^^ ~ 8.05159 and 6 is a constant. 

9.3. Parameters for connected SP-graphs. In order to approximate the growth 
constant for connected SP-graphs, we need to refine the analysis over the cycle in- 
dex sum for pointed 2-connected unlabelled SP-graphs, which is defined in terms 
of simpler pointed classes using the dissymmetry theorem for tree-decomposable 
structures developed in [S]. In particular, we have that Zqi ~ £' + Z, where 

Z — Z-jz + Zm ~ Z-jz^M ■ 

Here I' refers to a pointed link-graph (with cycle index sum si), and Z-jz, Zmi 
Z-R-M are the cycle index sums associated to the classes of pointed 2-connected 
graphs with a pointed 7?.-brick, a pointed A^-brick and an a pointed edge TZ — 
in the associated TZMT-tree (see [5] for proper definitions). Each one of these 
series can be written in terms of the cycle index sum associated to series networks. 
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general networks and networks which remains invariant when a change of the pole 
is applied. Denoting these series by Zg , Zxi and , respectively, we get 

(18) I Zm - ,si(2exp(E,>oiZ«) -2-2Z5-i(z| + zf )) 
\zn M = ■'*i^5(2exp(^,>oi4^)) -2-Z5) 

Recall that in the previous equation zl*-* = Z-^{si, S2i, ■ ■ ■)■ We denote by F{z) 
the solution of the equation F(z) = z e5cp{ZB'{F{z), F(z'^), . . .)) (exp denotes 
the Polya operator for sets in the unlabelled framework). This equation can be 
written in the following way: write J(z) — Zi3'{F{z), F(z'^), . . .) and A{z) = 

(Si>2 i"^ ('^')) ■ Consequently, F{z) = z exp ( J(z)) A(z). We define §(z) — 
Zs(F{z), F{z'^), ...) and E{z) = Zs{F{z), F(z'^), ...). Define also the following 
functions: E{z) = exp (e,>i M{z) = exp 27TT^ (^''^')) ' 

G{z) = exp ^X]i>2 (-^O) ■ Observe that all A{z), E{z), G(z) and M{z) are ana- 
lytic functions at the singularity of F{z) (which is the same singularity as the one 
of J{z)). With these definitions and using System (18) we get the following system 
of equations: 

J{z) = zexp(J(z))^(z) + 

(zexp(J(z))A(z) (2exp(§(z))G(z) - 1 - S(z)))^ (exp(§(z))G(z) - 1/2) + 
F(z^) (2exp(§(z2))G'(z2) _ 1 - g^^^)) (S(z)M(z) exp(§(z)) - 1/2) + 
zexp(J(z))A(z) (2 (exp(S(z))G(z) - 1 - §(z)) - {^zf + § (z^)) /2) - 
zexp(J(z))A(z)S(z) (2exp(S(z))G(z) - 2 - §(z)) 

S(z) = (2 exp(§(z))G(z) - §(z) - 1) z exp( J(z))v4(z) (2 exp(§(z))G(z) - 1) 

§(z) = (2exp(S (z2))G (z2) - 1) zexp(J(z))^(z) + 

(2exp(S {z^))G (z2) - 1) i^(z2) (2£;(z)M(z) exp(§(z)) - 1 - S(z)) . 

By a fixed-point argument we are able to compute from this system the first terms 
in the Taylor development of the series which appear in the previous equations. 
The procedure is the same way as in Section 9.2 (taking, for instance, the initial 
conditions A{z) = E{z) = G(z) = Af(z) = 1 and F{z) = §(z) = S(z) = 0). Later, 
we consider the simplified system in which we approximate all analytic functions 
by their Taylor series, up to a prescribed index. Consequently, the reduced system 
can be written in the compact form (J, §,§) = F(J, §,S;z). In order to find the 
critical points, we consider the associated singular system of equations, which is 
obtained from the system (J, S, S) = r( J, §, §; z) by considering the determinant of 
its Jacobian. 

Solving this system using a symbolic manipulator (for instance. Maple), we get 
that the singular value of z is p2 ~ 0.10655 for a truncation to order N = 50. 
Hence, the numbers Cn and of unlabelled connected and general SP-graphs on n 
vertices are 

where 72 « 9.38527 and c, g are constants. 
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